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SECTION  I 


INTRODUCTION 

The  macroscopic  properties  of  a  low  density  cloud  or  swarm  of  electrons 
in  a  neutral  gas  under  the  influence  of  an  electric  field  can  be  calculated 
from  solutions  of  the  Boltzmann  equation  or  with  Monte  Carlo  techniques  if  the 
appropriate  cross  sections  are  provided.  Information  on  such  properties  is 
critical  to  the  understanding  of  phenomena  in  gaseous  electronics.  When  the 
cross  sections  are  uncertain,  comparisons  of  calculated  swarm  parameters  with 
measurements  using  estimates  for  the  cross  sections  can  lead  to  refined  esti¬ 
mates  for  these  cross  sections.  Some  of  the  best  estimates  of  the  low  (ther¬ 
mal  to  several  volts)  energy  electron  scattering  cross  sections  from  atoms  and 
simple  molecules  have  been  derived  from  such  techniques.^  The  accuracy  of  the 
calculated  swarm  parameters  or  the  cross  sections  so  derived  depends  on  the 
accurate  solution  of  the  Boltzmann  equation. 

Many  approaches  have  been  used  for  the  mmerical  solution  of  the  Boltzmann 

2 

equation  for  equations  in  the  context  of  swarms.  By  far  the  most  common  meth¬ 
od  is  the  two-term  expansion  solution;  the  electron  velocity  distribution  func¬ 
tion  is  expanded  in  spherical  harmonics  and  the  expansion  is  truncated  after 

3 

the  first  two  terms.  Computer  codes  to  implement  this  technique  have  been 

4 

developed  and  widely  circulated,  especially  throughout  the  laser  community. 

Last  year,  under  this  contract,  a  multi-term  spherical  harmonic  expansion 
method  was  developed^  which  is  especially  suited  for  implementation  on  the 
large-scale  vector  computers  becoming  available  now.  Here  we  use  that  method 
to  Investigate  the  validity  of  the  two-term  approximation  in  N,.  We  have  also 
used  the  multi-term  method  to  study  the  effects  of  ionization  and  anisotropic 
scattering  cross  sections. 
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In  Section  II,  we  outline  the  multi-term  method.  A  detailed  description 
cun  he  found  in  the  final  report  from  1979.^  In  Section  III,  we  present 
studies  of  the  convergence  of  the  swarm  parameters  using  the  Boltzmann  code 
and  the  effect  of  boundaries  by  comparing  the  Boltzmann  results  with  Monte 
Carlo  studies  in  N,  at  moderate  field  strengths.  The  calculations  in  that 
section  were  done  using  a  set  of  cross  sections  previously  derived  from  swarm 
experiments'7  and  the  fields  were  constrained  to  bo  low  enough  that  secondary 
electrons  resulting  from  ionization  were  safely  ignored.  The  extension  of  the 
calculations  to  higher  fields  required  values  of  higher  energy  cross  sections 
than  were  available  in  the  swarm  derived  data  set  and  slight  modifications  of 
the  ca lcul at ional  techniques.  The  effect  of  anisotropic  scattering  was  inves¬ 
tigated  and  found  to  be  negligible  at  low  field  strengths  but  quite  important 
at  the  one  example  shown  for  higher  fields.  These  items  are  discussed  in 
Section  IV  and  we  present  calculations  up  to  field  strengths  high  enough  that 
about  ‘SO"  of  the  electron  energy  loss  Is  to  ionization. 


The  implications  of  this  work  are  discussed  in  Section  V. 


SECTION  II 


MATHEMATICAL  MODEL  AND  SOLUTION  METHOD 

The  Boltzmann  equation  for  electrons  defines  the  electron  energy  distri- 
♦  + 

bution  function,  f(r,v,t)  and  is  written  as, 


+  v*V  f  +  a*V  f  =  C(f)  (1) 

3t  r  v 

where  a  is  the  acceleration  due  to  an  external  electric  field  and  C  is  the 

collision  operator.  Since  Eq .  (1)  depends  on  three  variables,  r,  v,  and  t, 

some  assumptions  must  obviously  be  made  in  order  to  solve  for  f.  In  the 

2a  b 

regime  of  drift  tube  experiments,  the  usual  assimptions  ’  are 
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at 
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(2a) 
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(2b) 


These  two  assimptions  reduce  the  dependence  of  f  from  seven  variables  to  the 
three  velocity  variables.  At  this  point  a  further  assumption  is  usually  made 

■f 

namely,  that  the  angular  dependence  of  f(v)  can  be  approximated  by  the  first 
two  terms  of  a  spherical  harmonics  expansion,  or  rather  Legendre  expansion  be 
cause  of  the  cylindrical  symmetry, 

f ( v )  =  f( e ,  9)  =  fQ(e)  +  fj(t)  cos  0  .  (2c) 


The  first  assumption  [Eq.  (2a)]  is  met  satisfactorily  in  the  drift  tube  ex¬ 
periments  when  field  strengths  are  low  enough  that  there  is  no  appreciable 
ionization  and  where  there  is  no  attachment.  Then  the  electrons  do  reach  an 
equilibriun  fairly  quickly  on  the  time  scale  of  the  experiments.  The  second 
assumption  [Eq.  (2b)]  is  also  fairly  accurate  in  these  cases,  but  in  order  to 


calculate  diffusion  coefficients,  it  is  necessary  to  include  the  effect  of  the 

8  9 

spatial  gradients  in  some  way.  Parker  and  Lowke  and  Skullerud,  for  example, 
have  proposed  methods  for  doing  this.  The  third  assumption  [Eq.  (2c)]  is  the 
"two-term"  approximation.  In  order  for  the  two-term  assumption  to  be  valid, 
the  melastic  cross  sections  must  be  small  compared  to  the  elastic  cross  sec¬ 
tions,^1  i.e.,  on  the  average  the  electrons  must  lose  only  a  small  fraction  of 
their  initial  energy  upon  colliding  with  the  neutral  gas  atoms  or  molecules. 
This  is  not  always  the  case,  especially  for  molecules.  Implicit  in  the  two- 
term  approximation,  but  not  usually  stated,  is  the  idea  that  the  electron- 
neutral  cross  sections  are  no  more  anisotropic  than  cos  0.  In  other  words,  if 
the  cross  sections  were  also  expanded  in  spherical  harmonics,  it  would  only  be 
tiie  first  two  terms  that  would  enter  into  the  calculation. 

Tlie  method  we  use  to  solve  the  Boltzmann  equation  have  been  described  in 
detail  previously^’k  for  the  case  of  no  ionization.  We  will  summarize  the 
method  here  and  address  the  extension  to  the  case  where  ionization  is  present. 

■f 

We  include  the  t  and  r  dependence  of  the  distribution  following  the  method 
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proposed  by  Skullerud  in  which  the  full  distribution  function,  including  the 
r  and  t  dependence,  is  expanded  in  powers  of  the  spatial  gradient  of  the  elec¬ 
tron  density,  n(r,t); 


f  (  r ,  v,  t)  =  f  ^°\v)n(r  ,t)  -  f^\v)  •Vn(r,t) 

+  1  f^(v)  x  (-$)^n(r,t) 

k=2 


(3) 


The  first  term  in 
tion  function  and 


this  expansion  is  the  simple  product  of  a  velocity  distribu- 
a  density  distribution  function.  The  succeeding  terms  can 


be  thought  of  as  corrections  to  this  simple  product  approximation 


X 


Combining  this  expansion  with  the  Boltzmann  equation  we  write, 

f(o)(0)  ML&l  -  f(1)(0)  .  ^  mLi! 

v  '  at  v  '  at 

+  a*$vf  ^\v)n(r,t)  -  a  *$v  {?^  ^  (v)  *$n(r ,  t )  } 

+  f(0)v-tf  n(r,t)  -  v.tf  { v)  *^n(  r ,  t )  }  +  ... 
r  r 1  ' 

=  C{f  ^^(v)n(r  ,t)  -  ^  *$n(r  ,t)  >  +  ...  .  (4) 

The  continuity  equation  for  electrons  provides  an  expression  for  the  time 
derivative  in  the  first  two  terms, 

MLzI  =  w(0)n(r ,t)  -  ^(1)^n(?,t)  +  w(2)  :  ^n(?,t)  (5) 

3 1  “ 

where  the  u/^  is  the  electron  gain  frequency  due  to  ionization,  is  the 

(2) 

drift  velocity,  w  is  the  diffusion  tensor,  and  the  succeeding  w’s  are 
higher  order  transport  coefficients. 

Combining  Eqs.  (4)  and  (5), 


(0)£(0)  ,  e(0)  ^/£(0)Nl 


+  a  •  V  f 
v 


-  C(f  )  ]n 

+  [-u,(0)?(1)  +  a.^/(0  -  C(?(1)>]  •  $n  +  [v-w(1))f(0).^n 
+  . . ,  =  0 


(6) 


A  hierarchy  of  equations  results  from  setting  coefficients  of  powers  of 
the  gradient  of  n  equal  to  zero.  The  first  two  equations  are, 


a.vvf(0)(v)  -  J0)f(0)  -  C(f(0))  -  o 

a.v  r(1)(v)  -  J0)l(l)  -  C(?(1))  =  -(v-:(1))f 


(0) 


(7a) 

(7b) 


and  tne  u/s  for  k  >  0  can  be  written  as 


(8) 


(k)  ,  +  (k-l)3 

u>  =  J  vf  d  v  . 

These  equations  can  be  solved  in  order,  first  for  f^\  then  for  the  higher 
order  terms.  Thus,  the  effect  of  spatial  and  temporal  gradients  are  treated 
as  exactly  as  necessary  for  the  oj's  of  Interest.  The  only  iteration  required 
is  in  the  determination  of  in  the  first  equation  and  we  return  to  this 

point  later.  The  problem  is  now  reduced  to  solving  the  first  equation,  (7a). 
Hie  same  method  may  then  be  applied  to  the  second  and  succeeding  equations  if 
those  are  desired. 

(0)  + 

The  method  described  below  is  shown  for  f  (v),  but  the  same  methods 

+(1)  + 

have  been  applied  to  the  calculation  of  f  (v).  We  first  expand  the  distri¬ 
bution  function  in  spherical  harmonics,  or  rather  Legendre  functions,  because 
of  the  cylindrical  symmetry, 

N 


f(0)(v)  i  f(e,0)  =  i  fi(e)Pi(cos  6) 


(9) 


However,  rather  than  considering  only  two  terms,  we  retain  an  arbitrary  number 
of  terms  at  this  point.  The  substitution  of  the  spherical  harmonic  expansion 
into  the  Boltzmann  equation  results  in  a  set  of  coupled  differential  equations 
for  the  coefficients,  f^(v).  The  complication  in  solving  these  equations  comes 
from  the  nonlocal  nature  of  the  collision  terms  due  to  the  inelastic  and  super- 
elastic  collisions.  When  a  simple  two-term  expansion  is  used,  the  resulting 
set  of  coupled  differential  equations  may  be  combined  into  one  second-order 
differential  equation  for  fg(v)  which  is  often  solved  by  the  "backward  prolon¬ 
gation"*^  technique.  Rather  than  attempt  to  extend  the  backward  prolongation 
idea  to  the  coupled  system,  we  have  used  a  global  approach. 

In  order  to  implement  a  global  solution,  the  coefficients,  f^Ce),  are 
further  expanded  in  some  set  of  known  functions,  cubic  B-splines,  S(e).** 
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In  our  case. 


fi(e)  =  l  C  S  (e)  +  61 
j=l  J  -1 


(10) 


Since  the  splines  do  not  form  a  complete  set  because  we  do  not  consider  an 
infinite  nunber  of  them,  there  is  a  small  error  term,  6^,  in  the  expansion. 

The  error  can  be  forced  equal  to  zero  at  each  point  on  a  set  of  grid  points 
spanning  the  velocity  range  of  interest  by  requiring  < S ^ ( e ) I  6 =  0  for  i  = 
0,2,...^  and  j  =  1,2,...N  .  This  allows  us  to  write  the  set  of  coupled  dif¬ 
ferential  equations  in  the  form  of  a  system  of  linear  algebraic  equations  which 
can  be  expressed  as  a  matrix  equation,  MC  *  B,  where  the  M's  are  functions  of 
the  cross  sections,  E/N,  and  velocity.  The  elements  of  B  contain  the  boundary 
conditions.  This  equation  can  be  easily  solved  by  any  of  a  number  of  matrix 
inversion  methods  to  yield  the  desired  coefficients,  the  C's.  Given  the  C's, 
we  can  construct  the  fj(e)'s  and,  hence,  the  total  distribution  function. 

For  the  calculations  reported  here,  we  used  typically  100  splines  to 
cover  the  energy  range.  The  energy  range  was  from  zero  to  the  point  where  the 
distribution  function  had  dropped  five  to  seven  orders  of  magnitude  from  its 
peak.  The  grid  defining  the  splines  was  unevenly  spaced  for  maximum  flexi¬ 
bility  and  the  integrals  were  performed  with  range-splitting  Gauss -Legend re 
quadrature  techniques.  All  the  nonphysical  parameters  were  varied  to  ensure 
that  the  final  solution  depended  in  no  way  on  the  high  energy  cut-off  or  num¬ 
ber  of  splines. 
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SECTION  III 


MODERATE  FIELD  STRENGTHS 

In  this  section,  we  discuss  the  convergence  of  the  distribution  function 
and  transport  and  rate  coefficients  calculated  using  the  method  discussed  in 
Section  II  in  the  case  of  N2  over  a  range  of  E/N  from  1  to  200  Td  and  show  com¬ 
parisons  of  the  two-term,  multiterm  and  Monte  Carlo  results.  Our  expectation 
is  that  these  results  are  very  typical  of  other  molecular  gases  which  have 
similar  electron  scattering  cross  sections. 

In  this  range  of  E/N  in  N2  average  electron  energies  go  from  well  below 
the  onset  of  the  vibrational  thresholds,  through  the  vibrational  maximum,  and 
fall  slightly  short  of  the  onset  of  the  electronic  thresholds.  The  cross  sec¬ 
tion  set7  used  for  these  calculations  is  shown  in  Figure  1.  These  cross  sec¬ 
tions  are  isotropic.  The  momentum  transfer  cross  sections  shown  in  the  figure 
is  the  total  momentum  transfer,  the  sum  of  the  elastic  and  inelastic  momentum 
transfer  cross  sections.  The  rotational  excitation  cross  sections,  for  which 
an  example  is  shown  in  Figure  1,  were  replaced  by  a  single  level  cross  sec¬ 
tion.  The  vibrational  and  electronic  cross  sections  are  shown  in  the  figure 
as  sums  of  the  individual  level  excitations.  In  the  calculation  we  considered 
the  individual  level  cross  sections  and,  with  ionization,  a  total  of  23  in¬ 
elastic  cross  sections  were  included.  We  restrict  ourselves  here  to  the  range 
of  E/N  low  enough  that  ionization  may  be  treated  as  an  energy  loss  mechanism. 

At  the  highest  value  of  E/N  considered,  200  Td,  about  3%  of  the  electronic 
energy  lost  goes  into  ionization.  For  higher  E/N,  the  electrons  born  in  the 
ionization  events  begin  to  significantly  affect  the  distribution  and  must  be 
included  as  will  be  discussed. 
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ELECTRON  ENERGY-eV 


Figure  1 


Electron  collision  cross  section  for  N9  cross 

sections  used  in  the  calculations  of  Section  III. 
Excitation  cross  sections  are  shown  as  sums  over 
many  levels.  In  the  calculations,  all  of  the  23 
individual  cross  sections  were  used  rather  than 
the  sums. 


1.  Convergence  of  the  Distribution  Function  and  Transport  Parameters 

Since  the  computational  problem  is  the  determination  of  the  spherical 
harmonic  expansion  coefficients,  the  f^'s,  we  wiH  first  look  at  those.  The 
normalized  f^(e)'s  in  ^  at  an  E/N  of  100  Td  (where  the  average  energy  of  of 
the  electrons  in  the  swarm  is  2.2  eV)  are  shown  in  Figure  2  as  a  function  of 
energy.  The  results  shown  are  the  first  four  coefficients  in  a  six-term  ex¬ 
pansion.  The  isotropic  component,  f q ,  is  the  largest  at  all  energies.  The 
higher  order  coefficients  are  smaller  but  still  significant  in  comparison  to 
f^.  This  is  especially  true  in  the  2  eV  region  where  the  vibrational  cross 
sections  are  large.  The  considerable  structure  in  the  2  eV  region  and  again 
around  7-8  eV  seems  to  reflect  the  onset  of  important  inelastic  cross  sections 
at  those  energies.  As  will  be  seen  later,  the  f^  component  of  the  distribu¬ 
tion  function  calculated  from  only  two-terms  in  the  case  shown  may  be  quite 
different  from  the  multiterm  result. 

However,  only  the  first  few  coefficients  enter  into  the  calculation  of 
the  measurable  parameters  of  interest  in  most  swarm  applications  and  not  the 
total  distribution  function.  For  the  calculation  of  those  measurables,  it  is 
really  not  important  to  have  the  full  distribution  function.  The  important 
thing  here  is  to  determine  how  well  we  have  calculated  the  first  few  coeffi¬ 
cients.  Due  to  the  coupling  in  the  equations,  the  solution  obtained  for  f^, 
for  example,  will  depend  on  the  number  of  terms  in  the  expansion.  Figure  3 

shows  a  comparison  of  f^  calculated  from  a  two-term  approximation,  a  six-term 

1 3 

approximation,  and  a  Monte  Carlo  technique.  This  is  the  same  case  as  seen 
in  Figure  2,  ^  at  100  Td.  The  six-term  and  Monte  Carlo  values  of  f^  compare 
very  well  at  all  energies  except  those  around  the  origin.  Since  these  are  two 
completely  different  methods,  the  conclusion  drawn  from  this  is  that  f()  con¬ 
verges  very  quickly  as  a  function  close  to  the  correct  solution  even  though  the 
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Figure  2.  Electron  energy  distribution  for  N2»  The  first 
four  normalized  Legendre  coefficients  of  the 
distribution  function  in  N£  at  100  Td  in  a  six- 
term  calculation. 
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NORMALIZED  f 
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higher  order  coefficients  in  Figure  2  are  still  quite  large.  A  four-term  cal¬ 
culation  was  also  done  for  this  case  and  yielded  an  f^  very  close  to  the  six- 
term  value  indicating  independently  the  convergence  of  the  six-term  f^. 

The  convergence  of  the  transport  coefficients,  drift  velocity  and  dif¬ 
fusion  coeff icients ,  and  rate  coefficients  is  connected  very  closely  to  the 
convergence  of  f^  as  can  be  seen  in  Eqs.  (8)  and  (9).  Figure  4  shows  the  con¬ 
vergence  of  the  transport  coefficients,  drift  velocity  (v^),  transverse  and 
parallel  diffusion  (D,j,  and  Dp)  coefficients  and  the  A  £  electronic  excitation 
rate  in  case  of  nitrogen  at  100  Td  as  a  function  of  the  order  of  the  solution 
or  the  nunber  of  spherical  harmonic  components  used  in  the  calculation.  From 
the  figure  it  can  be  seen  that  the  two-term  values  of  the  drift  velocity  and 
the  transverse  diffusion  coefficient  are  higher  than  the  higher-order  calcula¬ 
tions.  Conversely,  the  values  of  the  A-state  excitation  rate,  illustrative  of 
electronic  excitation  in  general,  and  the  parallel  diffusion  coefficient  are 
lower  in  the  two-term  than  in  the  higher-order  calculations.  The  calculated 
values  are  seen  to  converge  as  the  order  of  the  calculation  is  increased. 

Beyond  four  terms,  there  is  very  little  change  in  the  values. 

Calculations  similar  to  the  one  presented  in  Figure  4  were  carried  out  over 
a  range  of  E/N  from  1  to  200  Td  and  the  difference  between  the  two-term  results 
and  the  six-term  results  are  shown  in  Figure  5  as  a  function  of  E/N.  The  con¬ 
vergence  of  the  parameters  shown  in  Figure  4  for  the  seven  values  of  E/N  inves¬ 
tigated  is  similar  to  the  100  Td  case  seen  in  Figure  3.  We  consider  here  the 
six-term  values  to  be  the  converged  results.  As  for  the  100  Td  case,  the  values 
of  the  drift  velocity  and  the  transverse  diffusion  coefficient  in  the  two-term 

calculation  are  higher  than  the  converged  result.  The  difference  increases  with 

3  3 

increasing  E/N  to  a  maximun  and  then  decreases.  Similarly  for  the  A  £,  C  JI,  and 
4 

W  A  electronic  excitation  rates,  the  difference  increases  to  a  point  and  then 


13 


%  Change  from  2-Term  Value 


Figure  4.  Convergence  of  transport  and  excitation  coeffi¬ 
cients  for  N2  at  E/N  =  icf15  V-cm2  =  100  Td.  The 
convergence  as  a  function  of  the  order  of  the 
solution  for  the  drift  velocity,  transverse  and 
parallel  diffusion  and  A  state  excitation  in  N2  at 
100  Td. 
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E/N  (Td) 


Figure  5.  Errors  due  to  2-term  solutions  of  Boltzmann  equa¬ 
tion  for  N2.  Error  introduced  by  the  two-term 
approximation  when  compared  to  the  six-term  solu¬ 
tions  for  drift  velocity,  transverse  diffusion, 
and  certain  electronic  excitation  rates  as 
functions  of  E/N. 
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ilrv  reuses ,  hut  Ih'  re  the  two-term  values  for  the  excitation  rates  are  lower  in 

general  than  the  converged  results.  The  maximum  difference  in  each  case  occurs 

at  70  Td.  At  this  value  of  E/N,  the  power  lost  by  the  electrons  due  to  elastic 

collisions  is  at  a  minimum  due  to  the  onset  of  the  vibrational  cross  sections, 

3a 

and  based  on  the  criteria  of  small  fractional  power  lost  per  collision,  we 

would  expect  the  two-term  app. ox imat ion  to  be  worse  in  this  region. 

Tlie  general  trends  in  extending  the  order  of  solution  from  two  to  six  terms 
are  an  elongation  in  energy  of  the  distribution  function  as  reflected  in  an  in¬ 
crease  in  the  electronic  excitation  rate  and  a  decrease  in  the  drift  velocity  and 

diffusion  coefficient.  These  trends  are  seen  in  both  Figures  3  and  A.  The  elon- 

1 A 

gated  tail  results  from  a  straggling  effect;  it  is  only  electrons  that  have 
avoided  collisions  that  can  acquire  a  relatively  high  energy  from  the  field. 

These  electrons  will  have  velocities  directed  almost  parallel  to  the  field  and  the 
"almost  spherical"  or  two-term  approximation  will  not  include  these  electrons. 

Attempts  to  parametrize  the  difference  between  the  two-term  calculation 
and  the  multi-term  calculation  have  been  only  partially  successful.  One  such 
attempt  is  shown  in  Figure  6,  the  fractional  error  in  the  two-term  transverse 
diffusion  coefficient  as  a  function  of  the  ratio  of  the  energy  exchange  colli¬ 
sion  frequency  to  the  momentum  transfer  collision  frequency.  This  ratio  is 
very  nearly  equal  to  the  ratio  of  the  drift  or  directed  energy  to  the  random 
energy.  When  this  number  is  large,  the  distribution  function  converges  very 
slowly  and  we  would  expect  the  two-term  approximation  to  introduce  significant 

error.  About  all  that  can  be  said  from  this  plot  is  that  the  error  always 

3f 

falls  above  a  straight  line  with  a  slope  of  unity.  Lin  et  al.  also  give 
criteria  for  the  error  in  the  various  coefficients  due  to  the  two-term  ap¬ 
proximation.  Their  equations  give  a  line  for  this  case  that  falls  below  ours 
bv  about  a  factor  of  two. 
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Momentum  Exchange  Freq. 
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Random  Energy 


Figure  6.  Error  in  transverse  diffusion  coefficient.  An 

attempt  to  parameterize  the  two-term  error  in  the 
transverse  diffusion  coefficient.  The  fractional 
error  is  plotted  as  a  function  of  the  ratio  of  the 
energy  exchange  frequency  to  the  momentum  exchange 
frequency  for  three  cases,  N2,  a  methane  model  and 
a  model  atom. 


Boundary  Effects 


In  order  to  compare  experimental  values  of  the  swarm  parameters  with  those 
calculated,  we  must  determine  how  well  the  mathematical  model  approximates  the 
experimental  situation.  Thus  far,  we  have  applied  the  Boltzmann  equation  to 
the  calculation  of  swarm  parameters  and  have  determined  that  the  approximations 
arc  valid  and  that  the  solutions  are  converged.  It  remains  to  show  that  the 
Boltzmann  approach  is  an  appropriate  model  of  the  experiments.  The  Boltzmann 
approach  taken  here  cannot  take  into  account  the  effect  of  the  physical  bounda¬ 
ries  present  in  experiments.  A  Monte  Carlo  simulation  was  done  to  access  the 
importance  of  the  boundaries  on  the  calculated  values  of  swarm  parameters. 

We  have  performed  Monte  Carlo  calculations  in  using  the  set  of  cross 

sections  presented  in  Figure  1.  The  simulation  conditions  of  a  one  cm  drift 

1 7  -3 

distance  and  a  density  of  10  cm  at  100  Td  were  chosen  to  approximate  the 
experiments  of  Tachibana,  Levron  and  Phelps^  and  Urosevic.1^  For  this  calcu¬ 
lation,  single  electrons  were  released  normal  to  the  cathode  with  a  2  eV  start 
energy.  The  electrons  were  allowed  to  drift  through  the  tube  and  were  absorbed 
a  the  cathode  and  anode.  Upon  absorption,  a  new  electron  was  released  with 
the  same  start  conditions.  The  simulation  continued  until  the  electrons  en¬ 
countered  10^  collisions  with  the  neutrals.  A  record  was  kept  after  each  col¬ 
lision  of  the  type  of  collision  and  the  position  of  the  electron  in  the  drift 
tube  when  the  collision  took  place.  From  this  record  we  have  plotted  the  num¬ 
ber  of  excitation  events  per  distance  interval  as  a  function  of  distance  as 

3 

seen  in  Figure  7  for  the  v  =  0+1  and  A  Z  excitations. 

The  effect  of  the  cathode  and  anode  can  be  seen  in  the  figures  as  a  depar¬ 
ture  from  the  equilibrium  excitation  rates.  One  would  expect  an  energy  equilib¬ 
rium  to  be  established  after  the  distance  required  for  the  drift  and  diffusion 
currents  to  become  equal  or  D/W  which  for  this  case  is  about  0.02  cm.  This  can 
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Figure  7.  Monte  Carlo  calculation  of  excitation  for  [ N 2 1  = 
1017  cm~3  at  E/N  »  10_15  V-cra2  =  100  Td.  Monte 
Carlo  calculations  of  a  number  of  excitation 
events  as  a  function  of  distance  in  a  drift  region 
for  the  v  =  0+1  vibrational  excitation  and  A  z 
electronic  excitation.  Excitation  rates  may  be 
determined  from  such  data  as  discussed  in  the 
text . 


1 
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bo  viewed  as  the  characteristic  distance  to  approach  equilibrium  or  the  1/e 

distance  to  equilibrium.  This  is  roughly  the  distance  seen  in  the  figures 

necessary  to  reach  the  equilibrium  rates.  The  v  =*  0*d  excitation  starts  out 

higher  than  the  equilibrium  because  the  release  energy  of  2  eV  is  equal  to  the 

energy  at  which  the  vibrational  excitation  cross  section  is  a  maximum.  Simi- 
3 

larly,  the  A  I  state  excitation  with  a  threshold  energy  of  6.17  eV  cannot  take 
place  until  some  electrons  have  gained  enough  energy  from  the  field  to  reach 
the  A  state  threshold.  These  trends  in  the  excitation  near  the  cathode  do,  of 
course,  depend  on  the  release  energy  chosen. 

The  anode  effects  are  independent  of  the  release  conditions  provided  an 
equilibrium  has  been  established  at  some  point  in  the  drift  region.  Here  we 
again  see  boundary  effects  on  a  scale  corresponding  to  the  D/W  energy  relaxa¬ 
tion  distance.  There  will  be  a  deficiency  of  low  energy  electrons  near  the 
anode  because  of  the  reduction  in  the  number  of  backscattered  electrons  due 
to  absorption  of  electrons  at  the  anode. 

The  excitation  rates  calculated  from  the  Boltzmann  code  and  those  from  the 
Monte  Carlo  simulation  are  compared  on  the  right  of  the  figures.  Hie  agreement 
i s  good  for  the  v  =  0+1  transition  and  only  fair  for  the  A  state  excitation. 
However,  when  excitation  events  near  the  boundaries  (0.1  cm  on  either  end)  are 
excluded,  the  agreement  between  Boltzmann  and  Monte  Carlo  rates  agree  to  within 
0.5%.  We  conclude  from  this  excellent  agreement  that  the  Boltzmann  treatment 
is  valid  in  real  experimental  situations  when  boundaries  are  present  and  that 
the  most  accurate  comparisons  between  the  Boltzmann  calculations  and  the  exci¬ 
tation  experiments  will  be  for  those  experiments  which  focus  on  the  center  of 
the  drift  region  rather  than  the  entire  drift  distance. 
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SECTION  IV 


HIGH  FIELD  STRENGTHS 

The  extension  of  the  method  and  the  N2  ‘calculations  presented,  thus  far, 
to  higher  values  of  E/N  requires  additional  effort  in  two  areas;  development 
of  the  calculational  method  to  include  the  new  electron  produced  in  the  ion¬ 
ization  events  and  the  extension  of  the  set  of  cross  section  data  to  higher 
energies  with  information  on  the  differential  scattering  cross  sections. 

Rather  than  attempt  to  modify  and  extend  the  swarm  derived  cross  sections 
used  in  Section  III  to  higher  energies,  it  was  decided  to  review  the  literature 
and  assemble  the  best  possible  set  of  cross  sections,  including  the  angular 
dependences  relevant  to  our  calculations.  Because  this  cross  section  set  has 
not  been  adjusted  to  yield  calculated  transport,  excitation  and  ionization 
coefficients  in  agreement  with  experiment,  this  set  should  not  be  regarded  as 
a  recommended  set  for  gas  discharge  laser,  etc.,  calculations.  This  exercise 
was  instructive  in  that  it  showed  that  while  a  large  part  of  the  necessary 
differential  cross  section  data  are  available  in  the  literature,  there  are 
some  potentially  important  omissions  and  discrepancies. 

1.  Differential  Scattering  Cross  Sections  in  N^ 

The  most  convenient  form  of  expressing  the  angular  dependence  of  the 
cross  sections  for  use  in  the  multi-term  Boltzmann  code  is  as  spherical  har¬ 
monic  components  of  the  cross  sections, 

Q1(e)  =  /  P1(cos  6)Q( 0 ,e)dil  , 

and,  rather  than  entering  the  Q^’s  f°r  aH  the  cross  sections  in  tabular  form, 
it  is  easier  to  express  the  Q^s  in  terras  of  the  Q0's,  the  total  cross  sections. 
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The  angular  dependences  of  the  cross  sections  for  elastic,  rotational, 


vibrational,  and  electronic  excitation  and  for  ionization  are  presented  in 

18  18 
Table  1.  The  magnitudes  of  the  cross  sections  were  based  on  a  previous 

19 

swarm  analysis  for  the  low  energy  elastic  and  rotational  cross  sections, 

20  21 

beam  experiments  for  the  higher  energy  elastic,  vibrational  and  elec- 

22  2  ^ 
tronic  cross  sections  and  on  theory  for  the  ionization  cross  sections. 


2.  Boltzmann  Treatment  of  Ionization 

For  the  higher  values  of  E/N,  the  average  electron  energy  is  such  that 
ionization  becomes  one  of  the  most  important,  if  not  the  dominant,  inelastic 
energy  loss  channel.  It  is  then  necessary  to  . nclude  the  new  electrons  pro¬ 
duced  in  the  ionization  event  in  the  calculation  of  the  electron  energy  distri¬ 
bution  functions.  The  mathematical  model  presented  in  Section  II  can  Include 
this  effect  and  the  appropriate  equations  are  Eqs.  (7)  with  a  non-zero 

a •  V v f ( v )  -  u,(0)f(v)  =  C[ f ( v) ]  (11) 

where  the  collision  term  must  be  modified  to  include  the  scattering-in  contri¬ 
bution  of  the  new  electron  as  is  given,  for  example,  by  Thomas. ^ 

It  is  instructive  at  this  point  to  write  the  Boltzmann  equation  explicitly. 
After  the  spherical  harmonic  expansion  and  subsequent  projection  of  equations, 


i 

!  i-1 


dfi-l(e)  i-1 


dfi+l ( e>  ,  i+2 


4.  N  1/2,  ,  .  _  _  N  _  *  ,  .  N  0  m  d  r  2.  ,  , 
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for  i  =  0, . . .N 
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RATIOS  OP  SPHERICAL  HARMONIC  COMPONENTS  OP  DIFFERENTIAL  SCATTERING  CROSS  SECTIONS  FOR  N.  USED  IN  BOLTZMANN  EQUATION 


where  QT  is  the  total  cross  section,  is  the  elastic  momentum  transfer  cross 

section,  is  the  ith  inelastic  cross  section,  is  the  threshold  energy  of 

the  kth  inelastic  process,  is  the  ionization  frequency,  and  I  is  defined 

later.  as  used  here  equals  in  Eq.  (5)  and  is  the  exponential  growth 

constant  of  the  total  number  of  electrons  fn  the  gap  in  a  pulsed  Tovmsend  type 

0  / 

experiment  as  discussed  by  Thomas.  The  first  term  on  the  right  represents 
electrons  scattered  out  of  an  energy  element  de  centered  about  e  due  to  all 
collisions.  The  second  term  on  the  right  represents  electrons  scattered  into 
ur  from  all  elements  de+  centered  about  e+  that  are  connected  to  de  by  inelas¬ 
tic  collisions. 

When  ionization  is  included,  the  scattering-out  term  does  not  change.  The 
contribution  of  ionization  to  the  scattering-out  term  is  independent  of  any 
electron  production.  However,  the  new  electron  will  contribute  a  scattering-in 
component.  We  assume  that  the  primary  and  secondary  electrons  share  the  excess 
energy  of  the  primary  over  the  ionization  threshold,  in  a  ratio  q:l-q  where 

0  <  q  <  1.  Then  in  order  for  the  primary  to  scatter  into  de  it  must  start  with 
an  energy  e+  =  (e/q)  +  e^  before  the  collision.  Similarly,  for  the  secondary 
to  scatter  into  de  it  must  have  originated  in  an  ionization  collision  between 
a  neutral  and  a  primary  electron  of  energy  e+  =  [e/(l-q)]  +  e^*  Since  the 
rate  of  events  leading  to  the  scattering-in  of  electrons  into  an  element  de 
centered  about  e  is  e+NQ(e+)f ( e+) ,  the  scattering-in  contribution  from  ioniza¬ 
tion  considering  the  secondary  electron  production  is, 


1  S  +  Ei)NfS+  eiMf  +  e!) 


+  TRT  (TmT  +  Gi )Nf *TI^qT  +  eiWo^T  +  £i) 


(13) 
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Hi  porously ,  q  is  a  function  of  £.  However,  the  work  of  Tagashira  and 
of  Hayashi26  suggests  that  setting  q  =  1/2  at  all  energies  is  a  good  approxi¬ 
mation.  We  have,  therefore,  used  q  =  1/2  in  all  the  following  calculations  to 
reduce  the  variation  of  parameter  space. 

3.  Numerical  Results  Considering  Ionization 

The  effect  of  ionization  is  not  significant  for  E/N  less  than  about  200  Td 
because  below  that  value,  less  than  3%  of  the  electron  energy  loss  is  due  to 
ionization.  We  carried  out  calculations  including  ionization  at  300,  500, 

1000  and  2000  Td  and  compare  to  calculations  where  ionization  is  treated  as 
an  energy  loss  only.  These  calculations  are  based  on  the  cross  sections  pre¬ 
sented  in  paragraph  1  of  Section  IV  but  using  isotropic  cross  sections  only  in 
order  to  isolate  the  effect  of  ionization. 

It  can  be  seen  in  Eqs.  (12)  and  (13)  that  the  inclusion  of  the  secondary 
electrons  produced  in  the  ionization  event  adds  a  term  dependent  on  v^/N  and 
changes  the  form  of  the  ionization  scattering-in  contribution  in  the  collision 
operator.  The  sign  of  the  vi/N  term  and  the  dependence  on  f  is  the  same  as 
that  of  the  scattering-out  contribution  in  the  collision  operator.  Therefore, 
when  this  term  is  included,  it  appears  that  the  scattering-out  has  increased. 
The  ionization  modified  scattering- in  must  balance  this  additional  scattering- 
out-like  effect.  The  overall  effect  is  a  shift  of  the  distribution  to  lower 
energies  because  the  scattering- in  electrons  come  in  at  lower  energies. 

The  solutions  of  Eqs.  (12)  were  carried  out  iteratively  because  v^/N  is  a 
function  of  the  distribution  function,  f.  It  vas  possible  to  find  a  consistent 
solution;  i.e.,  a  solution  where  the  input  value  and  the  calculated  value  of 

v^/N  were  equal,  by  straightforward  iteration.  However,  a  more  efficient  way 

1 9 

to  arrive  at  a  consistent  solution  was  to  calculate  the  energy  balance.  We 
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define  the  energy  balance  as  the  ratio  of  the  power  gained  by  the  electrons  due 
to  the  field,  FPG,  to  the  power  lost  in  collisions,  CPL.  In  a  consistent  solu¬ 
tion,  this  ratio  should  be  unity.  Departures  from  unity  are  due  to  an  incon¬ 
sistent  value  of  v^/N  and  it  is  possible  to  use  the  energy  balance  to  guide  the 
iteration.  We  write 

FPG  =  CPL  +  IPL 

where  IPL  is  the  ionization  power  lost  and  note  that  we  have  shown  that 

IPL  =  (<e>  +  e.)vi  , 


where  <e>  is  the  average  electron  energy.  The  energy  balance  is 


EB 


FPG _ FPG _ 

CPL+IPL  =  CPL  +  (<£>+£  )v± 


(14) 


While  it  is  not  possible  to  obtain  a  new  value  of  v1  from  this  relation,  Eq. 
(14)  does  provide  a  good  guide  to  the  selection  of  the  next  value.  If  EB 
is  less  than  (greater  than)  one,  vi  is  too  large  (small).  With  some  practice, 
the  number  of  iterations  can  be  reduced  to  two  or  three. 

Calculations  at  these  higher  E/N  values  are  done  with  three  approxima¬ 
tions,  two-term  without  ionization,  and  two-  and  six-terms  with  ionization. 
Based  on  the  convergence  studies  shown  in  Section  II,  we  assumed  that  the  cal¬ 
culations  at  these  E/N  had  converged  with  six  terms  in  the  spherical  harmonic 
expansion.  Several  quick  checks  were  made  at  500  and  1000  Td  which  confirmed 
this  assumption.  Values  of  the  drift  velocity,  W,  transverse  diffusion  coef¬ 
ficient,  Dt,  normalized  ionization  frequency,  v^N,  as  well  as  fractional  ion¬ 
ization  power  lost  are  shown  in  Table  2  for  four  values  of  E/N.  The  results 
are  much  as  expected.  We  see  by  comparison  of  the  two-term  solutions  with  and 
without  ionization  that  the  distribution  is  shifted  towards  lower  energies 
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with  the  Inclusion  of  Ionization.  The  effect  of  more  terms  In  the  expansion 
is  to  shift  the  distribution  function  towards  higher  energies  and  bring  the 
results  slightly  back  toward  the  two-term  values  without  ionization.  Very 
large  effects  appear  in  the  diffusion  coefficients  and  this  is  due  to  the  in¬ 
crease  in  the  spatial  gradients  when  ionization  is  present.  The  spatial  de¬ 
pendence  is  included  as  discussed  in  Section  II. 

We  see  then  that  the  inclusion  of  ionization  is  critical  to  an  accurate  cal¬ 
culation  of  swarm  parameters  at  high  E/N.  The  error  introduced  by  the  two-term 
approximation  is  dwarfed  in  comparison  except  in  the  case  of  diffusion.  There, 
the  two-term  approximation  is  more  severe  than  the  neglect  of  ionization. 

4.  Anisotropic  Scattering  Effects 

As  mentioned  above,  anisotropies  in  the  cross  sections  up  to  the  order  of 
the  Legendre  functions  retained  in  the  expansion  of  the  distribution  function 
may  be  included  in  the  multi-term  Boltanann  formulation.  That  this  is  the 
case  can  be  seen  from  the  collision  operator; 

NQT(e)ef1(e)  -  N  l  (e+e^f^e+e^  J  Qk< e+ek , 0)P .(cos  0)dfl 
k  x 

for  i  =  0, . . .  ,N  . 

The  integral  in  the  scattering- in  term  picks  out  only  the  ith  Legendre  compo¬ 
nent  of  the  cross  section.  For  example,  in  the  two-term  expansion  the  momentum 
transfer  cross  section  comes  up  because  it  contains  the  isotropic  and  the  cos  0 
components  of  the  differential  cross  sections,  the  same  angular  components  in- 
eluded  in  the  two-term  distribution  function.  In  the  multi-term  analysis,  a 
convenient  form  for  representing  the  angular  distribution  of  the  cross  sections 
l  s 

Q1  (e)  -  /  Q(e,0)P1(cos  0)dfl 
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We  first  calculated  the  distribution  function  and  the  associated  swarm  pa¬ 
rameters  from  1  to  200  Td  for  average  electron  energies  from  0.3  to  4.5  eV  with 
the  anisotropies  in  the  elastic  and  vibrational  cross  sections  given  in  this 
section.  The  neglect  of  anisotropies  in  the  electronic  excitation  cross  sections 
is  well  justified  for  E/N  <  100  Td  when  only  a  small  fraction  of  the  electron 
energy  goes  into  the  electronic  channels  but  probably  not  a  very  good  assumption 
at  200  Td.  We  find  that  for  these  cross  sections,  there  is  little  or  no  dif¬ 
ference  in  the  calculations  of  the  measurable  parameters  from  1  to  200  Td. 

In  order  to  compare  the  calculations  with  and  without  these  anisotropies, 
it  is  necessary  to  determine  what  to  keep  constant  between  the  calculations. 

We  have  kept  the  elastic  momentum  transfer  constant  and  in  so  doing,  the  f 
and  fj  equations  differ  only  by  virtue  of  cos  6  components  in  the  inelastic 
cross  sections.  Because  the  vibrational  anisotropies  do  not  include  a  cos  6 
component,  the  two-term  with  and  two-term  without  anisotropies  are  identical 
for  the  1-200  Td  calculations. 

For  higher  values  of  E/N  with  the  correspondingly  higher  average  electron 
energies,  the  anisotropies  in  the  cross  sections  become  more  pronounced  and 
potentially  more  important  in  the  calculations.  We  have  calculated  the  dis¬ 
tribution  function  and  the  swarm  parameters  at  500  Td  where  the  average  elec¬ 
tron  energy  is  about  8.8  eV  using  the  full  set  of  differential  cross  sections 
shown  in  Table  3. 

Table  3  shows  our  results  with  four  approximations,  the  two-term  without 
anisotropies  and  the  two-,  four-  and  six-term  with  the  anisotropies.  Ioniza¬ 
tion  is  included  in  all  four  cases.  For  these  calculations,  the  elastic  momen- 
tun  transfer  was  again  kept  constant.  The  only  difference  between  the  two-term 
with  and  the  two-term  without  anisotropies  is  then  due  to  a  different  inelastic 
momentun  transfer.  The  overall  effect  of  the  addition  of  anisotropies  is  to 
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Table  3 

THE  EFFECT  OF  ANISOTROPIC  SCATTERING  ON  THE  CALCULATED  VALUES  OF 


TRANSPORT 

AND  RATE 

COEFFICIENTS  AT 

500  Td 

Two-term 

Isotropic 

Two- term 
Anisotropic 

Four- term 
Anisotropic 

Six-term 

Anisotropic 

W  (xIQ7  cm  sec  *) 

3.35 

3.54 

3.32 

3.31 

D^.N  (xlO77  cm  ^  sec  ^ ) 

4.44 

3.40 

2.98 

3.01 

-9  3  -1 

A  state  (xlO  cm  sec  ) 

1.59 

1.59 

1.53 

1.53 

C  state  (xlO  ^  cm~*  sec  ^ ) 

1.75 

1.75 

1.67 

1.67 

-9  3  -I 

Ionization  (*10  cm  sec  ) 

1.20 

1.13 

1.16 

1.26 

<e>  (eV) 

8.84 

8.86 

8.66 

8.72 
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shift  the  distribution  towards  higher  energies.  This  is  to  be  expected  if  the 
anisotropies  are  prirr  rily  in  the  forward  direction  as  the  directed  velocities 
tend  to  be  preserved  in  the  anisotropic  case  and  destroyed  in  the  isotropic 
case.  The  effect  of  higher  order  expansion  coefficients  as  seen  in  the  com¬ 
parison  of  the  two-  and  six-term  results  with  anisotropies  is  differ*-"  from 
the  Isotropic  case  and  will  depend  on  the  particular  differential  ci  <-.  ec- 
tions  involved.  The  effect  of  anisotropies  in  Eqs.  (12)  is  to  add  scattering- 
in  terms  in  the  higher  order  equations.  Whereas,  before  the  only  source  of 
the  higher  order  coefficients,  the  f^'s,  was  due  t0  the  field,  we  now  have 
scattering-in  contributions  to  those  coefficients. 

The  preliminary  conclusion  drawn  from  these  examples  is  that  when  the 
electron  energies  are  such  that  an  appreciable  number  of  the  electron-neutral 
collisions  are  anisotropic,  the  effect  must  be  included,  and  the  number  of 
terms  in  the  expansion  must  be  high  enough  so  that  the  anisotropies  can  be 
well  represented  by  the  same  order  expansion  of  the  cross  sections.  A  more 
complete  investigation  of  the  effect  of  anisotropies  at  high  E/N  would  be  de¬ 
sirable. 
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SECTION  V 


CONCLUSIONS 

The  original  purpose  of  this  work  was  to  develop  a  method  of  solving  the 
Boltzmann  equation  suitable  for  the  extraction  of  cross  sections  from  swarm 
experiments  for  cases  where  the  usual  two-term  approximation  was  of  question¬ 
able  validity.  The  criteria  for  the  validity  of  the  two-term  approximation 

a-e  low  field  strengths  and  small  ratio  of  inelastic  to  elastic  cross  sec- 
la 

tions.  An  arbitrary  ratio  of  inelastic  to  elastic  cross  sections  can  be  ac¬ 
commodated  with  the  multi-term  method  and  the  high  field  cases  can  be  treated 
when  proper  account  is  made  of  electron  production  during  ionization  and  ani¬ 
sotropies  in  the  scattering  cross  sections.  The  extension  of  the  Boltzmann 
analysis  from  a  two-term  to  a  multi-term  spherical  harmonic  expansion  method 
was  accomplished  last  year  under  this  contract  and  a  detailed  descriptipn  may 
be  found  in  Reference  5.  A  summary  is  presented  in  Section  II. 

This  year  we  have  made  a  study  of  the  convergence  properties  of  the 
Boltzmann  solution  and  calculated  swarm  parameters  in  the  case  of  at  moder¬ 
ate  field  strengths  (1-200  Td) .  The  error  introduced  by  the  two-term  approxima¬ 
tion  is  small  (~1%)  for  drift  velocities,  larger  for  the  diffusion  coefficients 
(~5%)  and  even  larger  for  the  excitation  rates,  30X  in  the  worst  cases.  The 
maximum  error  occurred  at  an  E/N  of  70  Td,  where  the  electrons  are  "seeing" 
the  largest  ratio  of  inelastic  to  elastic  cross  sections.  Convergence  of  the 
transport  and  excitation  coefficients  with  a  number  of  terms  in  the  expansion 
to  within  a  few  percent  was  achieved  by  four  terms  in  all  cases.  These  errors 
are  significant  for  the  extraction  of  cross  sections  from  swarm  experiments  if 
those  cross  sections  are  to  be  compared,  for  example,  with  beam  data.  If, 
however,  two-term  derived  cross  sections  are  used  in  a  two-term  calculation 
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of  swarm  parameters,  these  errors  may  be  tolerable  in  applications  other  than 
cross  section  determinations. 

Effects  of  electron  production  in  ionization  and  anisotropic  electron 
scattering  must  be  considered  in  high  field  calculations.  The  division  be¬ 
tween  moderate  and  high  fields  here  is  based  on  the  amount  of  electron  energy 
lost  in  ionization.  E/N  values  greater  than  200  Td  in  ,  at  which  point  37. 
of  the  electron  energy  goes  into  ionization,  are  considered  high  fields.  The 
inclusion  of  electron  production  into  the  Boltzmann  formulation  is  straight¬ 
forward  but  the  resulting  equations  must  be  solved  iteratively.  With  experi¬ 
ence,  the  number  of  iterations  is  small.  The  results  of  Section  IV  show  that 
the  neglect  of  electron  production  is  more  severe  than  the  two-term  approxi¬ 
mation  at  high  E/N. 

Anisotropic  scattering  is  potentially  very  important  at  high  E/N  but 
negligible  for  low  E/N  with  the  anisotropies  given  in  Table  1.  Anisotropies 
in  the  cross  sections  introduce  scattering-in  sources  in  the  higher  order 
equations  for  the  distribution  function.  The  issue  of  convergence  then  has 
two  aspects,  convergence  of  the  spherical  harmonic  series  representation  of 
the  differential  cross  sections  as  well  as  of  the  distribution  function.  In 
contrast  to  the  low  E/N  calculations,  the  effect  of  anisotropies  in  at 
500  Td  is  quite  significant. 

Superelastic  effects  are  included  in  the  multi-term  Boltzmann  formulation, 
but  sufficient  data  have  not  been  compiled  at  this  time  for  cross  sections  from 
excited  vibrational  levels,  for  example,  to  make  an  investigation  of  the  ef¬ 
fects  worthwhile. 

In  summary,  the  calculations  of  electron  transport  and  rate  parameters 
using  the  multi-term  Boltzmann  method  can  be  made  as  accurate  as  required  for 
the  iterative  extraction  of  cross  sections  from  swarm  experiments  over  a  wide 
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range  of  field  strengths  corresponding  to  a  nearly  thermal  distribution  at  one 


end  and  a  distribution  markedly  affected  by  electron  production  in  ionization 
at  the  other. 
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